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Fabio Taddei, Michele Governale, and Rosario Fazio 
NEST-INFM & Scuola Normale Superiore, 1-56126 Pisa, Italy 
(Dated: February 2, 2008) 

Within the scattering approach, we develop a model for adiabatic quantum pumping in hybrid 
normal/superconductor systems where several superconducting leads are present. This is exploited 
to study Andreev-interference effects on adiabatically pumped charge in a 3-arm beam splitter 
attached to one normal and two superconducting leads with different phases of the order parameters. 
We derive expressions for the pumped charge through the normal lead for different parameters for 
the scattering region, and elucidate the effects due to Andreev interference. In contrast to what 
happens for voltage-driven transport, Andreev interference does not yield in general a pumped 
current which is a symmetric function of the superconducting-phase difference. 

PACS numbers: 73.23.-b, 74.45.+C 



Introduction. Pumping consists in the transport of par- 
ticles obtained, in absence of a transport voltage, by vary- 
ing in time some properties of a mesoscopic conductor. 
If the time scale for the variation of the scattering ma- 
trix describing the conductor is larger than the transport 
time, then the pumping is adiabatic and the number of 
particles transferred per period does not depend on the 
detailed time evolution of the scattering matrix but only 
on geometrical properties of the pumping cycle 0. 

Adiabatic pumping has attracted a vast interest, and 
different aspects of this phenomenon have been addressed 
U HL as, for example, the counting statis- 
tics of the pumped current, the generalization to multi- 
terminal geometries and the question of phase coher- 
ence. The idea of adiabatic pumping has been combined 
with other phenomena typical of mesoscopic physics, like 
spin-dependent transport (if! llll Il2j |. Kondo physics [l3l 
Il4j . Luttinger- liquid phvsics|15j. and Quantum Hall 
effect 0. So far, there have been only few investigations 
of adiabatic pumping in normal/superconductor hybrid 
structures. Zhou|l7| has considered the pumped current 
due to the time-modulation of the superconducting cor- 
relations induced in the normal region. Wang et al.\ufy 
have studied the combined effect of pumping and An- 
dreev reflection in a system with only one single-mode su- 
perconducting lead, finding up to a fourfold enhancement 
of the pumped current due to the interplay of Andreev 
and normal reflection. The generalization to a multi- 
mode superconducting lead was done by Blaauboer|l9|. 

In this paper we explore the physics of adiabatic pump- 
ing in the presence of several superconducting leads. In 
particular, we want to study Andreev interference in 
adiabatic pumping. Andreev interferometers have been 
intensively investi gate d in the past both in the diffu- 
sive limit [U El El Hi, and in the ballistic 
oneEEE3 (f° r an extended list of references see, for ex- 
ample, Ref. E3)- ■"■ n a standard Andreev interferometer 
(as those considered in Refs. |U El El El El E3) 
transport is driven by an applied voltage. In the present 
work, we study the problem of an Andreev interferometer 



when transport is induced by adiabatic pumping. 

The paper is organized as follows: we start by deriving 
a formula for the current pumped through a normal lead 
in the presence of several other normal and superconduct- 
ing leads; then we apply the formalism to a fork-shaped 
structure which exhibits Andreev interference. 

Formalism. We consider a system consisting of P 
normal and Q superconducting leads connected to a 
generic scattering region characterized by its scattering 
matrix S (matrices arc in boldface). The different su- 
perconductors are described by a constant pair potential 
A m = |A m | exp (?$,„), where m = 1, ...,Q labels the su- 
perconducting leads. We note that when all leads are in 
the normal state we can write the scattering matrix as 



R(e) T'(e) 
T(e) R'(e) 



(1) 



where R is a P x P matrix containing the scattering am- 
plitudes between the normal terminals; R' is a Q x Q 
matrix containing the scattering amplitudes between the 
superconducting terminals; T is a Q x P matrix describ- 
ing the scattering between the normal leads and the su- 
perconducting ones; and T' is P x Q matrix describing 
the scattering between the superconducting leads and the 
normal ones. The energy e is measured with respect to 
the Fermi energy. Writing the scattering matrix as in 
Eq. Q makes evident that the system is equivalent to 
one consisting only of a normal lead with P modes and 
a superconducting lead with Q modes, each mode in the 
superconducting lead having its own pair potential A m . 
We can now write the scattering matrix for the hybrid 
normal-superconductor system in Nambu space: 



S = 



n ee (e) 



(2) 



where the submatrices are obtained composing the ma- 
trix Sn with the scattering matrix Sns of a perfect 
NS multichannel interface |22j- The latter is a diago- 
nal matrix of Andreev reflection amplitudes which can 
be written, under the Andreev approximation, as Sns = 



2 



a e~ **, where a is a diagonal matrix whose elements 
are exp[— iArcos(e/|A m |)] and e~ l * is a diagonal matrix 
of elements exp[— i$ m ]. As a result: 

U cc {e) = R{e) + 

T'(e)ae i ^R'*{~e)ae- i ^M(e)T(e) (3a) 
TZ hc (e) = T , *(-e)ae- i *M(e)T(e) (3b) 

with 

M(e) = [1 - i?'(e)ae i * R* (-e)ae _i *] _1 . (4) 
In Eq. (0), 7?-oc is a P x P matrix of nor- 
I 



mal (Andreev) scattering amplitudes between the nor- 
mal leads. Note that we have used the particle-hole sym- 
metry, which yields 7Z-hh(e) = ^-cc(~ e ) an d ^eh( e ) = 

-Kc(-f)- 

By means of the scattering matrix Eq. © , operating 
along the same lines of Refs. |l9|, |30|, we can write the 
charge pumped through any of the normal leads: 

Q n = ~fdx 1 dx 2 rMXi,x 2 ), (5) 

/^Normal leads 

where X\, and X 2 are the two pumping fields, and 



n,u(Ai, x 2 



ck 



' de J m | 8X l 8X 2 



(e)) r 



d(K he *(e)) nJ d(n he (e)) 



dX l 



dX 2 



(6) 



with /(e) being the Fermi Dirac distribution. Although 
Eq. © is valid at finite temperature, in the rest of the 
paper we will restrict ourselves to zero temperature. 

Andreev interferometers. We, now, apply the formal- 
ism developed above to an example of an Andreev inter- 
ferometer. The most simple system which allows us to 
investigate pumping with different superconductors con- 
sists of a 3-arm beam splitter schematically shown in Fig. 
^ For the sake of simplicity we consider a symmetric 
beam splitter, whose scattering matrix (in the normal 
state) is 



-siyl - 2 7 y/j yfy 
Vl a b 

V7 b 



(7) 



a 



where a = l/2(s 2 + — 27), b — 

siy/1 — 27), with Si — ±1 and < 7 < 



l/2(-« 2 + 
1/2. Lead 
is normal metallic, 



1, on the left-hand-side of Fig. ^ 
while the other two leads, denoted by u and d, are super- 
conducting with order parameters, respectively, equal to 
|A u |exp(z$ u ) and | A d | exp (i$ d ). 

The parameters to be varied in time are the strengths, 
U n and £7 d , of two additional ^-barriers placed in the two 
arms on the right side of Fig. ^ When all leads are in the 
normal state, the total scattering matrix Sn is obtained 
by combining the scattering matrix Sb given in Eq. l|7| 
with the transmission and reflection amplitudes of the 
two (5-barriers, namely 

t = - (8) 



and 



1 + ifiu/d 



. Pu/d 

' 1 + iPu/d 



(9) 



where /3 u / d = U n /d/ (Jivy), vf being the Fermi velocity. 
We choose as pumping fields /? u /d, i-e. the strength of 
the (5-barriers. We consider the sinusoidal week pumping 
limit: f3 u — (3 U + S(3 U shi(ujt) and /3d = (3d + Sftd sm(u)t — 

First we start by studying the case when the normal 
lead is tunnel-coupled to the rest of the structure, ob- 
tained by setting 7 — * 0, si = — 1 and s 2 — 1 in Eq. JJJ). 
The analytical expression for the charge pumped through 
lead 1, in leading order in 7, reads 



^i7 2 4sin(A$)^ + ^- 2toc ° S(A$) 



(Pu + Pd) 



COS 



2 (¥) 



(10) 



|A u |exp (i(|) u ) 




|A d |exp (i<|> d ) 

FIG. 1: Schematic picture of the Andreev interferometer, 
consisting of a symmetric beam splitter with 5-barriers added 
on the two arms where the superconducting leads are con- 
nected. The strength of the 5-barriers can be varied in time, 
and are used as pumping fields. 



3 



being A<I> = $ u — $d the phase difference, and A = 
irSf3 u S(3d sin 5(f) the area of the pumping cycle in param- 
eter space. The 7 2 dependence of the charge in Eq. l(TU|) 
is the expected one for transport mediated by Andreev 
reflection. It is interesting to note that Qf is an odd 
function of the phase difference, and that no pumping 
occurs at zero phase (A<£> = 0). The pumped charge 
Qf can be contrasted with the linear DC conductance 



Gdc = x 2 l C^he)i,i I 2 of the system 28], with the bar- 
riers frozen at their average values /3 u ,d, when a transport 
voltage is applied between the normal and the supercon- 
ducting terminals (superconductors being at the same 
potential). For this particular case, in leading order in 7, 
Gdc reads 



G 



DC 



e 2 2 1 + 4 (ft + ft) + 8 (ft + ft) + (1 + 4ft)(l + Aft) cos(A$) 



(Pu + Pd) 



(ii) 



The DC conductance, Eq. iJH}, is an even function of the 
phase difference. It has a zero-phase extremum, which 
can be either a maximum or minimum depending on 
the strength of the barriers. To complete the analysis, 
we mention that the system acts as a pump also when 
all leads are in the normal state. To leading order the 
pumped charge is linear in 7 (as expected for vanishing 
superconducting gap), and it reads 



Qi 1 



-A—4ry- 



Pu~Pd 



(ft + Pd 



1 



(12) 



In contrast to the superconducting case, the leading order 
of the pumped charge vanishes when ft = Pa . 

Now, let us turn to the case of a maximally- 
transmitting beam splitter, which is obtained from 
Eq. J7J) setting 7 = 1/2 and si_2 = 1- The analytical form 
for the charge pumped through lead 1 is rather involved 
and we report only the limits P u /d *C 1 (large barrier 
transmission) and P u /d 3> 1 (small barrier transmission) : 



Qf 



256Ai- (/3u 



Id) cos 2 (4/) 
"tt [3 + cos(A$)] 3 

sin(A<I>) 

n^ftpf 



if P u/d < 1 



if p u/d > 1 



(13) 

It is interesting to note that while in the case of large 
barrier transmission Qf is an even function of the phase 
difference, for small transmission Qf is odd. However no 
definite parity is present for arbitrary transmissions P n / d 
and, in particular, when /3d = we obtain: 



s _ 128Ai [1 



Qf = 



(1 + API) cos A<j) - p u sin Ac/)} 



[3 + 10pl + cos A0] 3 



(14) 



Again, it is instructive to see what happens for the DC 
conductance, when the barriers are frozen to their aver- 
age value. Also for this case, we report the same two 



limiting cases shown for Qf 
cos 2 (4/) 



G 



e 2 

X 32 



DC 



[3 + cos(A$)] z 

e 2 ft + ft + 2ftftcos(A<S>) 
h W4 



if ft /d < 1 



if Pu/d > 1 



(15) 

The DC conductance is an even function of the phase 
difference both large and small barrier transmission. It, 
actually, remains even also for arbitrary values of /3 u /dj 
while its zero-phase extremum can be either a maximum 
or a minimum depending on the values of /3 u /d- Finally, 
we report the pumped charge when the system is in the 
normal state 



-A-(p u 

TT 



Pd) 



if P n/d <l 



-A e - 



1 



if 

Pu Pd, 



(16) 



Pu/d > 1 



'nspipi 

Contrasting Eq. (|13fl with Eq. (|16(l . we notice that for 
the case /3 u /d < 1 , to leading order in the pumping 
parameters, superconductivity produces an enhancement 
of the charge pumped, reaching a maximum of a factor 
4. 

Finally, we wish to point out that the most distinc- 
tive signature of Andreev interference in the adiabatic 
pumping regime is the lack of a definite symmetry of the 
pumped current under inversion of the superconducting- 
phase difference. On the contrary, the DC current pro- 
duced by an applied transport voltage, either DC or time- 
dependent, is always an even function of A<F The case 
of a DC voltage has been discussed above. It can be eas- 
ily seen that also the current produced by rectification of 
an oscillating voltage V osc (for example induced by the 
pumping voltages on stray capacitances |3l|) is an even 
function of A<F In fact, the current produced by recti- 
fication reads ^J, |3a ^icct 



LU/2TT J^ /UJ G(r)V osc (r)dT, 



4 



where G(r) is the instantaneous conductance which is 
an even function of A$, and all other quantities do not 
depend on the superconducting-phase difference. This 
lack of symmetry with respect to superconducting-phase 
difference can be exploited to distinguish between pump- 
ing and rectification. This is analogous to the normal 
case where the symmetry used for this purpose is the one 
related to magnetic field inversion [sij . 

Conclusions. In this paper we have derived a scat- 
tering formula for adiabatically pumped charge in hy- 
brid NS multi-terminal systems. This has been used 
to study Andreev interference in a 3-arm beam splitter 
attached to one normal and two superconducting leads 
with different phases of the order parameters. Within 
the weak pumping limit we found that Andreev interfer- 
ence very much affects the charge pumped through the 
normal lead, though differently with respect to the case 
of DC-voltage-driven transport. In general, the pumped 
charge has no definite symmetry under inversion of the 
superconducting-phase difference and no zero-phase ex- 
tremum is found. 
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